Collective spontaneous emission from a system of quantum dots 
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We study the spontaneous emission from a regular lateral array or a randomly distributed ensemble 
of quantum dots under strong excitation (full inversion) conditions. We focus on the similarities and 
differences between the cases of random and regular arrangement of the dots and show that there 
is very little difference between the evolution of luminescence in these two cases, both for identical 
dots and for a realistically inhomogeneously broadened ensemble. This means that the enhanced 
emission or superradiance effect is not due to accidental clustering of pairs of dots. Moreover, we 
point out that observation of an enhanced emission under weak excitation does not prove that true 
superradiance will develop in a fully inverted system. 
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I. INTRODUCTION 

The phenomenon of collectively enhanced emission, de- 
scribed by Dicke for ensembles of atoms can also 
occur in quantum dots (QDs) 0. The cooperative ra- 
diation effect is due to the collective interaction of the 
two emitters with quantum radiation field. In the case 
of non-identical dots and in the absence of interactions, 
the collective emission appears if the interband transition 
energies in various QDs differ by no more than the emis- 
sion line width which requires the dots to be nearly 
identical, beyond the recent technological possibilities. 
However, coupling between the dots restores the collec- 
tive nature of the emission and leads to accelerated or 
slowed down emission even for dots with different transi- 
tions energies [J]. 

Two typical couplings that may appear in a system 
of QDs are long-range Coulomb interactions and short- 
range couplings that may result from a combination of 
carrier tunneling and Coulomb correlation effects. As we 
have shown, these two couplings show no qualitative dif- 
ference in either accelerating or slowing down the emis- 
sion in spite of their essentially different physical origin 
and properties [j| . The simulations of spontaneous emis- 
sion from ensembles of coupled QDs show that the emis- 
sion rate in such a system can indeed be increased due to 
collective coupling of the emitters to the electromagnetic 
field @. As the positions of the dots are random within 
the sample plane, it might happen that this enhanced 
emission effect results from clustering of dots in a self- 
assembled sample, as a result of which strongly coupled 
pairs of dots are formed and contribute to the emission. 
Therefore, in order to fully understand the role of the 
inter-dot coupling it is interesting to study the role of the 
planar arrangement of the dots on the observed collective 
emission. Moreover, so far the collective emission from 
QD ensembles was discussed only under weak excitation 
conditions (one excitation delocalized over the ensemble) 
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FIG. 1: The planar arrangement of the QDs: (a) regular ar- 
ray; (b) random ensemble (sample dimensions in nanometers). 



[g , where it is manifested only in the form of an enhanced 
emission rate. In order to observe actual superradiance, 
that is, a non-monotonic evolution of the emission with 
a delayed outburst of radiation Q , one has to study the 
photon emission under strong excitation conditions. 

In this paper, we study the collective spontaneous 
emission from small ensembles of QDs under strong exci- 
tation (full inversion) conditions for two kinds of spatial 
arrangements of the QDs in the sample plane: regular 
arrays and randomly scattered dots. We show that the 
dynamics of emission very weakly depends on the QD 
arrangement. 

The paper is organized as follows. In Sec. HH we de- 
scribe the model of a small ensemble of QDs in either ran- 
dom or regular planar arrangement. Next, in Sec. IIIIl we 
present and discuss the simulation results for the photon 
emission in these two cases. Finally, Sec. IIVI concludes 
the paper. 



II. MODEL 

We consider an ensemble of a few (up to six) QDs 
placed in the xy plane at positions r Q , where a num- 
bers the dots. Each QD is modeled as a two-level sys- 
tem (empty dot and one exciton) with the fundamental 
transition energy E a = E + e a , where E is the aver- 
age transition energy in the ensemble. We consider two 
different arrangements of the QDs in the sample plane: 
Regular arrays of QDs (double dots, triangular, rectan- 
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gular aligned structure of the dots, penta- and hexagonal 
structures, see Fig.QJa)) and ensembles of randomly dis- 
tributed QDs, as schematically depicted in Fig. QJb). In 
the latter case, the ensemble is randomly distributed with 
a fixed surface density over square mesas with the restric- 
tion that the inter-dot distance cannot be smaller than 
10 nm. The spectral properties of the dots are modeled 
by a Gaussian distribution of their transition energies 
E a with the standard deviation a around the mean E. 
The dots are coupled by an interaction V a p which can 
be either of dipole-dipole character (long-range disper- 
sion force) or short-range (exponentially decaying with 
the distance). 

We introduce the transition operators for the dots: the 
"cxciton annihilation" operators a a which annihilate an 
exciton in the dot a, and the '"exciton creation" opera- 
tors o\ which creates an exciton in the dot a. The ex- 
citon number operator for the dot a is then n a = a^ a a a . 
Using these operators, the Hamiltonian of the many dots 
system is written in the frame rotating with the frequency 
E/H in the form 

The long-range dipole coupling is described by 

r(lr) 



V a0 = V$ = -hT G(k r aP ), 
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is the spontaneous emission (radiative recombination) 
rate for a single dot, Eq is the vacuum permittivity, e r 
is the relative dielectric constant of the semiconductor, 
k = nE/ (he), where c is the speed of light and n = ^fel 
is the refractive index of the semiconductor, and 
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where f a p = r a p/r a p and d = d/d, where d is the in- 
terband matrix element of the dipole moment operator 
which is assumed identical for both dots. For a heavy 
hole exciton, d = (d /^/2)[l, ±i, 0] T , so that for a vec- 
tor r a f3 in the xy plane one has \d ■ r Q ^| 2 = 1/2. The 
short-range coupling is described by 



V 



(sr) 



V e 



The effect of the coupling to the radiation field is ac- 
counted for by including the dissipative term in the evolu- 
tion equations, which describes radiative recombination 
of excitons. The equation of evolution of the density ma- 
trix is then given by Q 
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FIG. 2: The interference term of the decay rate Fi2 and the 
short- and long-range coupling amplitudes VI r , V S r as a func- 
tion of the inter-dot distance. In (a), the small distance sec- 
tion is shown, while in (b) the oscillating tail at larger dis- 
tances is visible. 



where T aa = Fpp = T , T a p = T p a = T F(k r al3 ), 
with 
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and {...,.. .}+ denotes the anticommutator. The diago- 
nal decay rates T aa describe the emission properties from 
a single dot, while the off-diagonal terms T a p, a ^ (3, ac- 
count for the interference of emission amplitudes result- 
ing from the interaction with a common reservoir and are 
responsible for the collective effects in the emission. The 
values of the two couplings as well as the interference 
term of the decay rate T a /3 are plotted as a function of 
the distance between the dots in Fig. [2j 

In our simulations, we use the parameters for a typical 
CdSe/ZnSe QD system: T = 2.56 ns" 1 , n = 2.6, the av- 
erage transition energy of the QD ensemble E = 2.59 eV 
and the QD surface density v = 10 11 QDs/cm 2 . For the 
tunnel coupling we choose the amplitude Vq = 5 meV 
and the range r = 15 nm. 



III. RESULTS 

In Fig. [3] we show the results of numerical simulation 
based on Eq. [T] for regular QD arrays. At the initial 
time, the system is strongly excited to the fully inverted 
state \ip(0)) = |l..l..l) and we subsequently determine 
the photon emission rate as a function of time. On each 
plot, the photon emission in the system is shown as a 
function of time for identical dots (a = 0) and for non- 
identical dots with the realistic value of the fundamental 
transition energy standard deviation a = 18.4 meV 

As can be seen in Figs. GHa) and (b), the sponta- 
neous emission from a regular array of identical QDs de- 
velops a maximum, which is a few-emitter counterpart 
of the delayed outburst of photon emission observed in 
macroscopic samples Q- This effect appears already for 
three QDs and becomes more and more pronounced as 
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FIG. 3: The evolution of photon emission for the fully inverted 
initial state in regular arrays of a few QDs. 
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FIG. 4: The evolution of photon emission for the fully inverted 
initial state in an ensemble of a few QDs with a random spatial 
distribution. 



the number of dots increases. It implies that the true 
superradiance effect appears in the regularly arranged 
small ensemble of identical QDs with either long- or 
short-range couplings. However, if the array is formed 
by non-identical dots with a realistic inhomogeneity of 
the fundamental transition energies (Fig. [UJc) and (d)), 
when the value of a is greater than both the coupling 
strengths and the inverse exciton life time (the emission 
line width), then the maximum is completely washed out. 
Thus, under conditions that led to enhanced emission in 



the weak excitation regime [6|, full superradiance effects 
under strong excitation are not observed. 

In order to find out whether a randomly distributed 
ensemble of QDs behaves similar to a regular array of 
QDs we now study the photon emission from a random 
QD ensemble in the same strong excitation regime. In 
each case, we calculate an average of 100 realizations of 
the evolution in systems with a given number of dots but 
differing in their positions within the sample plane. Figs. 
EJa) and (b) show the simulated photon emission for an 
ensemble of identical dots. Again, a superradiant maxi- 
mum builds up as the number of QDs grows. Similarly 
to the case of regular QD arrays, the maximum disap- 
pears and the decay of luminescence becomes monotonic 
when the realistic distribution of fundamental transition 
energies is taken into account (Figs. UJc) and (d)). 

By comparing Fig.[3]with Fig.0]one can see that there 
is very little difference in arranging the dots either ran- 
domly or regularly under strong excitation conditions. 
Regular ordering of QDs only slightly enhances the peak 
emission in the hypothetical case of identical dots and 
slightly suppresses the emission in the realistic case of an 
inhomogeneously broadened ensemble. 



IV. CONCLUSIONS 

We have shown that the delayed outburst of radia- 
tion, typical for the superradiant emission develops in 
the luminescence from regular arrays or randomly dis- 
tributed ensembles of quantum dots in the strong excita- 
tion regime. The way the dots are distributed in the sam- 
ple plane (regular vs. random) makes very little differ- 
ence on the photon emission. This means that the system 
response is not dominated by accidental clustering that 
might appear in the random distribution case and lead 
to strongly enhanced contribution to the overall emis- 
sion from pairs of accidentally very closely spaced dots. 
Moreover, in both cases, the superradiant maximum is 
completely washed out if the realistic degree of inhomo- 
geneity of the fundamental transition energies is taken 
into account. Thus, enhanced emission observed experi- 
mentally under weak excitation does not imply that true 
superradiance will be manifested for a fully inverted sys- 
tem. 
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